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Electroweak radiative corrections are computed for Higgs production through vector boson fusion, 
qq — >■ qqH , which is one of the most promising channels for detecting and studying the Higgs boson 
at the LHC. Using soft-collinear effective theory, we obtain numerical results for the resummed 
logarithmic contributions to the hadronic cross section at next-to-leading logarithmic order. We 
compare our results to HAWK and find good agreement below 2 TeV where the logarithms do not 
dominate. The SCET method is at its best in the high LHC energy domain where the corrections 
are found to be slightly larger than predicted by HAWK and by other one-loop fixed order approx- 
imations. This is one of the first tests of this formalism at the level of a hadronic cross section, 
and demonstrates the viability of obtaining electroweak corrections for generic processes without 
the need for difficult electroweak loop calculations. 
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I. INTRODUCTION 

One of the major goals of the Large Hadron Collider 
(LHC), is the discovery of the Higgs boson and the study 
of its properties. The dominant production mechanism 
is Higgs production through gluon fusion, gg ^ H . Vec- 
tor boson fusion (VBF), qq — > qqH, is the second largest 
production channel and its cross section is about an or- 
der of magnitude smaller. Nevertheless, its measurement 
is very important for constraining the Higgs couplings 
and thus identifying the nature of the Higgs sector, see 
e.g. Refs. [l|,i. 

In this paper we will determine the electroweak cor- 
rections to the VBF process, using soft-collinear effec- 
tive theory (SCET) The framework 0, Q for 
calculating electroweak radiative corrections in SCET 
allows one to obtain the logarithmic contributions 
a" 2 In™ Sij/M^ ^ to the cross section, where Sij are the 
(partonic) generalized Mandelstam variables. When 
are parametrically of the same size and Sij 3> M^^, 
these logarithms dominate the electroweak corrections 
and may even require resummation [ol-fTTj. We will 
obtain numerical results at next-to-leading logarithmic 
(NLL) order, using the extension of the above frame- 
work to VBF in Ref. [l^l, and including the effect of 
parton distribution functions (PDFs). This is one of the 
first tests of this formalism at the level of a hadronic 
cross section. In our knowledge only a previous test has 
been reportedfisj. quite recently, regarding the top-quark 
forward-backward asymmetry in the process qq — >■ tt. 

The tree- level diagrams for VBF are shown in Fig.l. 
The two outgoing quarks produce forward jets with a 
large rapidity gap, which characterizes the VBF chan- 
nel. This allows one to suppress reducible and irreducible 
(i.e. gg ^ H + 2 jets) backgrounds using selection cri- 
teria called VBF cuts. These consist of tagging the two 
forward jets and sometimes a veto on central jets, see 



e.g. Refs. [TM3. The VBF cuts restrict us to a region of 
phase space that does not satisfy the above assumptions 
on the Sij too well. We explore the validity of our results 
in this region by comparing with HAWK [3 [13 ; and find 
good agreement below 2 TeV. At the higher energies of 
LHC the method is expected to be more reliable, because 
the Mandelstam variables are quite larger than the elec- 
troweak scale and the large logarithms dominate the elec- 
troweak corrections requiring some resummation which 
is here provided by SCET through the renormalization 
group (RG) runnning of the coefficients in the effective 
interaction. A comparison with HAWK shows that in the 
7-14 TeV range of LHC the electroweak corrections are 
larger than predicted by fixed order perturbative approx- 
imations. A subtle problem arises in this energy range, 
in some phase space regions where the Mandelstam vari- 
ables happen to be quite different in size, because the 
choice of the high energy scale seems to be ambiguous. 
The same problem also arises in one-loop and tree-level 
calculations. Our choice of an average scale seems to in- 
terpolate smoothly between phase space regions where 
the logs are small and regions where they dominate the 
corrections. 

This case study suggests that Refs. [1, US] may be used 
to directly obtain electroweak corrections for generic pro- 




FIG. 1: The tree-level t, u and s-channel diagrams for elec- 
troweak Higgs production. The s-channel process is usually 
referred to as Higgs strahlung rather than vector boson fusion. 
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cesses, even in regions of phase-space where they are not 
exphcitly vaUd, without the need to perforin difficult elec- 
troweak loop calculations that would have to be carried 
out for each individual process. These electroweak cor- 
rections are obtained in analytical form, and can be easily 
inserted in the software packages that have been devel- 
oped for computing QCD corrections to cross sections. 
We will discuss in some detail how to combine these elec- 
troweak corrections with known QCD corrections, which 
is essentially multiplicative at low orders in perturbation 
theory. 

QCD corrections to the VBF process have been cal- 
culated at next-to-lea ding order (NLO) for the inclusive 
cross section in R ef. [211 and for the differential cross 
section in Refs. p^l22ll23j. The gluon- induced contribu- 
tion at next-to-next-to-leading order (NNLO) was deter- 
mined in Ref. [IJl and the structure-function approach 
was employed in Refs. \2E-^ to get an accurate approx- 
imation to the full NNLO QCD corrections. The NLO 
electroweak corrections were determined in Refs. [l8l[l9l|. 
and are comparable in size to the QCD corrections and 
thus numerically important. 

We will now discuss in more detail how electroweak 
corrections are obtained in the framework of Ref. 0, Q . 
Here we only consider the standard model gauge group. 
Extensions like the minimal left-right symmetric gauge 
group f28l - l30| will be the subject of an other paper. The 
first step consists of matching onto SCET at a high 
scale fill ^ \/ Sij ■ This matching can be done in the 
unbroken SU{5) x SU(2) x U{1) gauge theory, since 
symmetry breaking effects are suppressed by powers of 



v/y^Jij. Using the renormalization group evolution, one 
then runs the effective theory operators down to a low 
scale /i; ^ Mz- At the low-scale, the W and Z boson 
are integrated out and one matches onto a SU{3) x U{1) 
effective theory, which only contains gluons and photons. 
The effects of SU{2) x U{1) symmetry breaking only en- 
ter in this low-scale matching. 

For VBF the amplitude is explicitly proportional to 
the vacuum expectation value (VEV), so the effective 
field theory operator is not a gauge singlet and standard 
resummation methods do not apply. Thus VBF provides 
a very interesting test of the method. The extension 
to the VBF process was derived in Ref. and the 

expressions obtained there will be used through out this 
paper. 

The paper is organized as follows: in Sec. [IT] the de- 
tails of the calculation are given, starting with the kine- 
matics in Sec. Ill Al followed by the high-scale match- 
ing in Sec. IIIBl the running in Sec. Ill C[ the low-scale 
matching in Sec. Ill Dl and the electroweak cross section 
in Sec. Ill El The necessary Passarino-Veltman functions 
are given in App. |^ and the Higgs wave function and 
tadpole contribution are given in App. [BJ In Sec. Ill Fl 
we discuss how to combine electroweak corrections with 
known QCD corrections, but postpone the correspond- 
ing numerics to future work. Our numerical results are 
shown and discussed in Sec. IIIII We conclude in Sec. IIVI 



II. CALCULATION 



A. Kinematics 



We will start by discussing the kinematics for the VBF 
process, qq qqH. The momenta of the incoming 
quarks (beams) are denoted by pi and p2 and the mo- 
menta of the outgoing quarks (jets) are denoted by and 
P4. The Higgs boson momentum is ph and we will assume 
that it is produced on-shell, = M^. We will neglect 
quark masses, pi = p| = pi = pI = 0, and follow the 
conventions of Ref. [l2l | by taking all momenta to be in- 
coming. The kinematic configuration is fully determined 
by the following generalized Mandelstam variables: 

s = (Pi + P2f , h = (pi + Pzf , ti = {P2 + Pif , 

U3^{p2+P3f, U4 = {pi+P4f. (1) 

The other combinations can be related to these 

s' = (P3 + Pi)^ = Ml - S - <3 - i4 - "3 - "4 , 

(pi + Phf = Ml + 2pi -ph = Ml - s - h - Ui , 



(P2 + Phf = Ml - S - <4 - U3 . 

(pa + Phf = s + t4 + u4, 

(P4 +Phf = S + <3 + U3 . 



(2) 



In SCET, each colhncar direction has a set of two light- 
cone reference vectors rii = ±(l,ni), hi = ±(1,— n^) 
associated with it, where in our conventions we take 
the plus (minus) sign for incoming (outgoing) parti- 
cles. For the quarks we can simply take rii = ^Pi/Pi- 
For the Higgs we have to take its mass into account, 
Uh — — P/i/V (Phf - Ml- The product Ui ■ nj is then 



= ± 



Pi -Pj 

pIp' 



ni-n/i = ± 1 - 



Pi-Ph 



PiPh 



Ml) 



(3) 



with the plus sign when the particles are both incoming 
or both outgoing. The pi ■ pj and pi ■ ph can be directly 
written in terms of Eq. ([T]), using also Eq. ([2]). 

At hadron colliders, such as the Tevatron and LHC, 
the colliding quarks (i = 1, 2) move along the beam axis 
and carry a fraction Xi of the hadron momenta. 



Pt 



T E — 



(1,0,0,±1) 



(4) 



where E^m is the center-of-mass energy. In the cross sec- 
tion, the momentum fractions Xi are integrated over. The 
distribution of momentum fractions carried by the quarks 
in the proton is described by PDFs. We parametrize the 
momenta of the outgoing quarks by 

P3 = -£'3(l,sin6'3,0,cos6'3) , 

P4 = —£'4(1, sin 6*4 cos if, sin 6*4 sin cosOi) , (5) 
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where Ei > is energy and 9i the angle with the beam 
axis of particle i. The 9i are related to (pseudo)rapidities 
rii by tan(0i/2) — exp(— ry^). The azimuthal angle be- 
tween the particles is given by ip. Momentum conserva- 
tion fixes the momentum of the Higgs and yields one ad- 
ditional constraint, as is clear from the three-body phase 
space, 



d$a 



^3^^ (2^)32i?. (2^) 
X {2nfS'^{pi +P2+ P3 +P4+ Ph) 

/ n [E,AE,A{cosei)\d^e 



X ,5 
1 



i=3,4 
4 



Mi 



d(cos 6'3) d(cos O4) dip dE-- 



/I 



(6) 



In the first step we carry out the ph integral using the mo- 
mentum conserving delta function, and write pi in terms 
of the spherical coordinates in Eq. ([5]) . In the second step 
we use the on-shell condition for the Higgs to perform the 
E4, integral, yielding £'4 = /1//2 with 

fl = XiX2E^^^^ ~ Ml ~ [Xi +X2 + {X2 - Xi) COS03]EcmE3 



/2 = [Xl + X2 + {X2 - Xi) COS 94]Ec 

— COS ip sin 03 sin 04]i?3 . 



2\1 



COS 63, COS 04 

(7) 



The remaining integrals in Eq. ([6]) will be carried out 
numerically. Their boundary conditions are 



< 03,04 < TT, 



Ml 

EL 



<xi<l, 



< < 27r , 

< a;2 < 1 , 



Ml 



0< E-,< 



X1X2EI 



Ml 



[xi + X2 + {X2 — Xl) COs6'3]i?c: 



(8) 



where we included the bounds on momentum fractions of 
the PDFs. The boundaries will be further restricted by 
any cuts we impose on the final state. 

We can express the above variables in terms of the 
Mandelstam variables in Eq. ([T|), 



Xl 



Er 



X2 



1 . 

^3 =-:rT=(*3e 



E, 



E-i — - 



-Y 



cm 



2^' 



cos ^3 



smip 



-U3 e 

Y 

Use'' 



E4 



2Vi 



t4 — M4 e 
t4 A- U4 e"^ 
s' 



U4 e 

-Y 



1 



cos 6*4 



sm f 3 sm f 4 \ 
where we also needed 



Use' 

[ 1 — cos 0s cos 04 



2EsE4 



(9) 



Y 



1 Xx 

o log — 

2 X2 



(10) 



The rapidity Y describes the boost of the partonic center- 
of-mass, which can of course not be expressed in terms 
of the Lorentz invariant Mandelstam variables. 
For future reference, we include 



n,.p, = ±2pO, n,,-pn^~pl^^{plY~Ml, (11) 
with a plus (minus) sign for incoming (outgoing) quarks. 

B. High-scale matching 

1. Operator basis 

At the high scale, we can ignore the effect of elec- 
troweak symmetry breaking and work in the unbroken 
phase of SU{2) x C/(l), since the partonic center-of-mass 
energy y/s is large compared to the VEV. The basis of 
SCET operators for vector boson fusion is given by [l^ 

OiA.B = OiOa,OiOb, 
OsA,B = OsOa , OsOb , 

04A,B.c = OlO\ , 010% , OtO^ . (12) 
The Higgs sector is described by Oi , . . . , O4 



'2 

Os - $100 - , 

01 = $Ir"(/.o + cl)lT-<fh , (13) 

where $h — W^^(f>n^ denotes the coUinear scalar doublet 
(j)n^ that will produce the Higgs, plus the corresponding 
coUinear Wilson line W„^. The field (/jq denotes a soft 
scalar that will attain a VEV in the broken phase. The 
basis of operators for the quarks is 



Oa = *37''r"«'i*47M?^"*2 , 

Ob = Cf'I'37^*i*47m*2 , 

o^ = *37''r»vi-i^4^^vi>2, 

0% = *37''*1*47a.T"*2 , 



(14) 



where the subscript i — 1 , . . . , 4 on the field labels the 
particle with momentum pi. All these fermion fields con- 
tain coUinear Wilson lines ^'i = W^.ipm, as required by 
coUinear gauge invariance. In Eq. (|14l) . is a fermion 
doublet (singlet) if it is left-handed (right-handed). We 
suppressed the projectors Pr^l = (1 ± 75)/2 to keep our 
notation general and allow for both left- and right-handed 
quarks. We can consider each helicity separately and 
combine the contributions at the end. Operators Oa and 
are only well-defined if all quarks are left-handed, and 
O'X and 0% are only allowed if at least one of the quarks 
is left-handed. 
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matching 


running 




high scale 


low scale 


7 r /? 


LO 


0-loop 




1-loop 


NLO 


1-loop 




2-loop 


LL 


0-loop 


0-loop 


1-loop 1-loop 


NLL 


0-loop 


1-loop 


1-loop 2-loop 2-loop 


NNLL 


1-loop 


2-loop 


2-loop 3-loop 3-loop 



TABLE 1: Order counting in fixed-order and resummed per- 
turbation theory. F and 7 denote the cusp and non-cusp part 
of the anomalous dimension, and P denotes the beta function. 



The other Wilson coefficients vanish at tree level. Here, 
the variable ry^, is defined as 77^ = 1 if the particle with 
label i is left-handed and rji = if the particle is right- 
handed. The hypercharge Y — l/Q for left-doublets, 
y = 2/3 for right-handed up-type quarks and Y ~ — 1/3 
for right-handed down- type quarks. The contribution for 
the u-channel and s-channel can be obtained through a 
permutation, as discussed in Sec. Ill El 



C. Running 



The basis of operators in Eq. (ITi|) is not complete, be- 
cause it assumes that the incoming particles are quarks 
and because it only suffices for the t-channel contribu- 
tion. If one (or both) of the incoming particles is an anti- 
quark, we can simply obtain the corresponding basis and 
expressions by interchanging 1 o 3 or (and) 2 o 4. The 
PDFs will of course differ, and the PDF for an anti-quark 
is suppressed compared to the quark case. Similarly, to 
switch to the u-channel we interchange 3^4. We can 
thus obtain the results for the other channels by sim- 
ply making replacements and will therefore restrict the 
discussion to the t-channel with incoming quarks. Only 
when we square amplitudes to obtain the cross section 
in Sec. Ill El will we need to be careful in combining the 
different channels. 



The running of the Wilson coefficients is described by 
the RG equation 



(17) 



where 7 is the anomalous dimension for the operators. 
It should be noted that this is a matrix equation, corre- 
sponding to the 10 operators in Eq. (fT2|) . In Ref. p^ . 
the anomalous dimension was derived, 



7 + 7 
a 
An 



T,,-To(4Lh + 2)+^T,-To 4L, 



(18) 



2. Tree-level matching 

We now perform the high-scale matching, which can 
be done in the unbroken phase of the electroweak gauge 
theory. Since we are working up to NLL order in the 
electroweak corrections, the tree-level matching suffices, 
as shown in table HI The relevant terms in the Lagrangian 
that couple the scalar doublet to the gauge fields are 



The full-theory diagrams were shown in Fig.l. Depend- 
ing on the parton types, only some of the diagrams con- 
tribute. Phenomenologically, vector boson fusion is most 
interesting when the jets are in the forward direction, 
for which the s-channel contribution is suppressed. We 
therefore do not include the s-channel. Matching the t- 
channcl diagram for incoming quarks onto the operators 
in Eq. yields 

n I \ *3|(m) r' ( \ 2i.9]'(M) 



2^3*4 

CiAW = — —7 mY2 , 

CiBW - —7 -»1»72 ■ 

13 14 



(16) 



The 7 denotes the anomalous dimension one would 
obtain by summing over the soft and collinear func- 
tions [1, [l^l for all the external particles (i.e. the quarks 
and Higgs, but not the VEV), and 7 describes a new con- 
tribution. In Eq. (|18p. the sum on i runs over the quarks 
and we introduce the shorthand 



L,; = l0£ 



m - Pi 



(19) 



(15) The gauge group generator acts on particle i as, 



(20) 



The anomalous dimension was derived in the unbroken 
phase of the gauge theory, hence the presence of the gen- 
erator To that acts on the soft field (po (in the broken 
phase 00 gets a VEV). 

We now give the explicit expressions, needed for our 
numerical analysis. Under SU{2), the ten operators in 
Eq. (dSl) split into subsets {Oia, Cis, C2A, C2B, C2c} 
and {O^A,O^BiOiA,OiB,Oic} ^ that do not mix under 
renormalization. In these bases, the SU{2) anomalous 
dimension is identical for both subsets and will be given 
in terms of 5 x 5 matrices. The SU{2) soft function is 



5 



SU{2) 
75 



r(a2) 



3 , 



/ -4ci + 2c3 

C2 

C2- 

V2(C1+ - C2+) 



3C2 



3ci_ 

3C2- 





C2- Cl_ 
Ci_ C2- 

-4?7i[/i3 + 2ci+ C2 

C2 -4r72?724 + 2c2^ 



-2(c3_-ci_) 



Cl+ - C2+ 



C3+ - C2- 
-4ci + C3 + Ci+ + C2+/ 



(21) 



where 



constants are given by 



— rt, -71., — iO 
log^^^^ 



(22) 



and 



Ul2 
U23 



C/34), 
C/34), 



C2 = rjimiUii + U23 

C3 iiim{Ui2 + Ui4 

ci± = JJilC^ih ± Ush) , 

C2± = ??2(C^2/i ± Uih) , 

C3± = vimiiUu " c/34) ± (t/23 - U14)] 



(23) 



Note that if only one of the rji is zero, Eq. (1211) contains 
entries for operators that are not ahowed when one of 
the fields is an SU{2) singlet, e.g. Ci in the upper-left 
corner. However, they are harmless since the Wilson co- 
efficients for these operators vanish and these entries do 
not lead to mixing with any of the allowed operators. 
The cusp anomalous dimension V{a) fsil is known up to 
three loops [s^. We only need the two- loop result 



r(a) 



1+E 



f- 

V47r 



K 



in) 



V9 3/ 9 9 



(24) 



We obtained the scalar contribution from the fermion 
case, since both only enter through the vacuum polariza- 
tion at this order. For SU(2), np and ns are the number 
of fermion and scalar doublets. For U{1)y, npTp and 
nsTs get replaced by the sum over the squared hyper- 
charges Y of the particles. Explicitly, the group theory 



SU{2) :CA = 2,TF = Ts = -,nF = 6,ns = l, 
Uil)Y:CA = 0,nFTF^5,nsTs^l. (25) 



We now move on to the U{1)y soft function, which 
mixes the Higgs operators in Eq. but not the quark 
operators in Eq. (|14l) . This mixing takes place in the sub- 
sets {01,03} and {O2, O4} of Eq. p3)) . and the anoma- 
lous dimension for each subset is identical. In this basis 
the U{1)y soft function is given by 



U{1) 

75 



r(«i) 




where 



C4 = Y^^Uu + YiU24 + YiY2{Uu + U23 - Ui 

C5 = YiiUih - Ush) + Y2{U2h - U4h) ■ 



(26) 



C/34), 
(27) 



The collinear functions for the quarks and for the Higgs 
are given by 



7c 



lb 



-r(a2)??» + r(ai)F/ 



9^2 



|r(a2) + ^Tia,) 



3q!2 Cti 
: — 



3ai..2 



-Y^ 



3y? 



An 



47r 167r2 



(28) 



where we included the contribution from the top Yukawa 
yt to the Higgs wave function renormalization. 
The new contribution in Eq. (fT8|) is given by [l3| 



6 



7 



.SU{2) 



7 



T{a2)Lh 



Q2" 

2^. 



T{ai)Lh 



fx 



27r. 



/-3 




Vo 

l2x2 



o\ 

-3 
10 
10 
l/ 




?72i2/4 

3771L1/3 -2?7iii3 

3/72-f'2/4 

2772L2/4 



-ir(ai)(yiLi/3 + 1-2^2/4) 



»7l-^l/3 

V2L2/4 



-2772^24 

2?7i-^3/i 



»yi '72-^13/24 ^ 



'72^2/4 
V1L3/I 
-771772^1234/ 



(29) 



where we have used the same bases as for the soft anoma- 
lous dimensions and abbreviated 

Lab.../cd... ^ {La + Lb + ..:)-{Lc + Ld + ...). (30) 

The last line in Eq. ([29|) is the Higgs rescattering contri- 
bution, with coefficient = {3, 0, 0, 1} for the operator 

{0l,02^ 03,02}- 

The total anomalous dimension in Eq. (|18p is then 
7 = 7f +7s^^^^ +E 7c +7^^^^' +7^^^^- +7^ , (31) 



where each of the ingredients has to be appropriately 
written as a 10 x 10 matrix and the sum on i runs 
over the quarks and the Higgs. A simple crosscheck on 
this expression comes from reparametrization invariance 
(RPI) [Hill]. Under the RPI-III transformation, 

< ^ e'^'Tif , nf ^ e-'^-nf , (32) 

where Hi is a parameter that can be chosen different 
for each coUinear sector. The individual ingredients in 
Eq. (|5T|) are not RPI-III invariant but we find that their 
sum is invariant, as should be the case. 



D. Low-scale matching 

1. Operator basis 

At low energies the effects of electroweak symmetry 
breaking need to be taken into account. In this section we 
will match onto a basis of operators in the broken phase 
of the gauge group. The Higgs part of the operators gets 
matched onto a v{pLi)hn, where we include the running 
of the VEV. The quark part in Eq. (|14l) can be matched 



onto the set 

Oa = U37^'"1"47m"2 , 

Ob = U3^^uid4jf^d2 , 
Oc = d3l^diU4'yf^U2 , 
Od = d3j'^did4jf_,d2 , 
Oe = d3j''uiU4jfj,d2 , 
Of = W37^di J47^U2 . (33) 

Here u and d denote up and down-type fields. For op- 
erators Oa,---,0_d, each pair of fields 'tpj^tl; is a pair 
of left-handed or right-handed fields, whereas in Oe and 
Op all fields are left handed. 



2. Tree-level matching 

We first perform the tree-level matching for the Higgs 
part of the operator, given in Eq. At tree level, the 

soft scalar field simply attains a VEV and the collinear 
scalar field produces a Higgs, leading to 

{Oi, 0§, O3, 02} ^ {1, 0, 0, ~S''y2}vhr, . (34) 

We now match {Oia,b, 04a,b,c} onto the quark opera- 
tors {Oa, . . . , 0_f} in Eq. ([55]) . ignoring O2 and O3 which 
vanish at tree level. This is described by a 6 x 5 matching 
matrix 

O, ^ Vhr. Rf , C.^vY, Rf C^, (35) 

J i 

such that 

^0,0, ^^C./i,,©,. (36) 

i i 
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The matrix i?*^"^ is given by 

/ 1 3Cu,lCu,2 
-1 3Cu,lCd,2 
-1 3Cd,lCu,2 
3Cd,lQ,2 



1 

2 
V 2 



~Cu,2 
-Cd,2 
C«,2 
Q,2 






Cu,l 
-CdA 
Cd,l 










-1 

1/ 



(37) 



Some of the quark fields in the unbroken basis in Eq. (jl4p 
had to be left-handed, but some can be either a left- 
handed doublet or a right-handed quark of a specific fla- 
vor. To cover all these scenarios we introduce, 



Cu,i = 77i -I- (1 - rii) {Vi + 
Cd.i = ??i + (1 - 'ni){^ - Yi) 

3. One-loop matching 



(38) 



In the one-loop corrections to the matching, we need to 
take into account that only the W and Z boson (but not 
the photon) are integrated out, and we need to include 
the mass-splitting between the W and the Z. (Since 
the W and Z mass only enter in the one-loop low scale 
matching, we can neglect their running.) The following 
relation will be useful, 

a2T,-T, -l-aiY.Y, (39) 
= avi/(T+T7 + TrT+) + azTf Tf + acmQ^Qj- 
= a2(T,-T, - TfT^) + (a2 + ai)Tf Tf + acmQ^Qj , 

where T± = (T^ ± iT2)/^/2, = - sin^ OwQ, Y is 
the hypercharge and Q is the charge operator. 

In Ref. the low-scale matching was written as 

a pure quark contribution plus additional contributions 
that involve the Higgs sector. We first discuss this quark 
contribution, which is given in terms of the soft and 
collinear functions of Ref. [1, [l^l, and which we write 
below as operators acting on quark fields. This contribu- 
tion vanishes for the antisymmetric Higgs operators Of 
and O3 in Eq. ([T3l) . so we can use the same basis as in 
the tree-level matching, {Oia.s, C4A,b,c}- 

Using Eq. ([M)) . the soft function for the quarks in the 
low-scale matching is given by ^ 



Ds^^ log ^ E(-TrT, + TlT^n, 



2^ 



(40) 



where the sum is over pairs [if) of quarks. We can eval- 
uate the Ti • Tj part of Dg in the unbroken phase of the 
gauge group, 

SU{2) 

D's = -y] TrT,U,j = with ci±, C2± ^ , 

i-[0i2) 

(41) 



(ij) 



and then convert it to the broken basis in Eq. p3p using 
i?*^"-' . The remaining terms in Eq. (PO)) can only be eval- 
uated in the broken phase, so we first multiply by i?'*^' 
to match onto the broken basis. In the broken basis, 

(y) 

= -diag(ci-fC2,Ci-C2,Ci-C2,Ci-|-C2,C3-Ci,C3-Ci) , 

iv) 

= diag(5^ iC/13 + gus9u,2Ce + 91,2^24 

gl,lUl3 + guA9d.,2CQ + 5d,2C^24 
9X1U1Z + 9d,igu,2C6 + gl^2U24 

gl.iUiz + gd..igd,2CQ + gj9U24 



glUi4 + gugdcj + gjU23, 
glU23 + gugdcj + glUu) , 



(42) 



where gu,i and gd^i denote the coupling of quark i to the 
Z boson, 



1 2.2 

gu,i = -77i--sm Ow , gd,'. 



i?7, + i sin^ 6*^ . (43) 



(For Oe and Op all the quark fields are left handed, so 
we suppressed the subscript "i".) In Eq. we used 
some of the variables Ci from Eq. (|23| and introduced 



(44) 



C7 = Uiz + U24 - U12 ~ U34 . 



iij) 



The collinear functions for the quarks in the low-scale 
matching are 

^ ^^D{Mw,ni-p^,^i) + '^gfD{Mz,n^-p.,,^i) , 

(45) 

which is evaluated in the broken basis. Here, gi is gu,i 
{gd,i) for an up-type (down-type) quark and 



^ We did not separate the cusp and noncusp anomalous dimen- 
sion in the low-scale matching, since we only need the low-scale 
matching at one-loop order (see table 



D M, A* = 2 log - log — - - log2 — 
^ /i^ 2 

3, AP 57r2 9 
- o log ^ - ^ + 7 • 
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(46) 
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The pure quark contribution to the low-scale matching 
at one loop is thus 



Dq = vhn 



M, 



^logi^(i?(°).i?^-i?^i?(°)) (47) 

2 



giogi;;fi?f.i?('')+exp(5:i?^)i?(°) 



In the formulation of Refs. [8|, |20|, the full low-scale 
matching is exponentiated. However, only the collinear 
functions contain large logarithms, so the soft function 
does not need to be exponentiated to sum the large log- 
arithms. This is a convenient simplification because only 
the soft function contains off-diagonal matrix elements. 

We now give explicit expressions for the additional con- 
tributions to the low-scale matching derived in Ref. [l^] . 
For OiA,Bi hi addition to the quark contribution de- 
scribed above, we get 



Di = vhr, 



+ 



A 



Mlv 



IGtt' 



O, 
(48) 



where O = Oa,b is the quark part of the operator. This 
contribution is simply an overall factor and requires a 
factor of i?'"' to convert to the broken basis. The finite 
part of the Higgs wave function i?'' and the finite contri- 
bution from the tadpole T*^ are given in App. [BJ The re- 
maining contribution comes from rescattering, where the 
sum on J runs over {h, W, Z} with combinatorial weight 
ryj = {3,2, 1} and Bq is the finite part of the Passarino- 
Veltman fimction defined in Eq. (jAll) . 

For 02A.B,c we do not have a pure quark contribution, 
due to the antisymmetric combination of scalar fields. 
The only contribution is 



Do 



vh„ 



"^lo.^^(t/,.-flog^)(nO" 



47r 



Ml 



-T?03)-f^logi;^^(c/,:. + log^)TfO^ 



(49) 



where 0° = 0\ b c the quark part of the operator. 
We can evaluate the T^O^-piece in the unbroken phase. 



D' 



E 



Uih + log 



TfO° 




mLl/3 C2- 



mLi/3 ci, 

■ ^72-^2/4 



C2+ + ??l'72^13/24 





(50) 



which is a 2 X 3 matrix in the basis {Oa,b} x {0'\ g q}, 
and can be converted to the broken basis using _R(°) . The 
ci± and C2± were defined in Eq. (j23l) . The remaining 
terms in Eq. (|49l) are evaluated in the broken phase. First 
we convert —0\ g (j/2 to the broken phase by multiply- 
ing with [Note that the factor of -1/2 for 0% in 
Eq. ([M)) was absorbed in R^^\] In the broken basis. 



log^H? 



= \ diagfyyics -I- 772C9, rjiCs - 772C9, 



E 



- ?7iC8 -I- ?72C9, -?7iC8 - ?72C9, 

mm{cw - cii), -rjir]2{cio ~ cn)) 
(L/., + log^)Tf 



diag(g„,iC8 -I- g„,2C9, 5„,iC8 -I- gd,2Cg, 

gd,lCs + 5u,2C9, gd.lCs + gd,2C9, 

guCw + gdCii, gdCio + guCu, (51) 



where 



C8 = Uih — Ush + Li/3 , 
ClO = Ulh — Uih + il/4 : 



Cg — U2h — Uih + -^2/4 



Cll 



U2h - U- 



3/i 



^2/3 ■ 



(52) 



It is worth pointing out that Eq. (|49| is actually /i- 
independent, as is clear in Eqs. (f50| and (|5T|) . 

For O^A.B we do not have a pure quark contribution 
either. The only contribution is 



L'3 



, az. M. 
vh„ — lo, 
47r 



^'^T.iu.H + log^yfO. (53) 



This can be evaluated by converting the quark part of 
the operator to the broken phase using i?*^"^ and then 
applying Eq. ([5T|) . 

For OiA,B,c, in addition to the quark contribution in 
Eq. gZl), we get 



Di 



vh 



+ h \o 



A 



J2vtBo{M,,Ml) 



U,h)i€''''^VO'' 
(54) 
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The combinatorial weights for the rescattering contribu- 
tion are now 77^ — {3, —2, 1}. The new contribution is 
given by 



We evahiate the last term in Eq. ((54|) in the unbroken 
phase, 



h = 



aw 
An 



(2L,, + l)log 



M, 



w 



- 1 



(55) 



where 



/s(«;,l) = l + ^-2W^tanh-\/^4 (56) 
3 V w V w — 4 



2Li2 



;+V(i^-4) 



ui ui— ■\/ (w — 4)u; 
2Li2 



Ci+ - 7/iLi3 (-C2- + ?/2i2/4)/2 

C2+ - 772^24 (Ci_ + 7;iL3/i)/2 

-C2-+772L2/4 Cl-+»/li3/l (Ci+ + C2+ - 7?lLl3 - 772^24)72^ 



(57) 



This matrix is given in the basis {0\ b c) ^ b ci- 
We can then convert ~0\ b cl^ broken phase by 

multiplying with R^^^ . As for Di, D2 and I?3, the explicit 
/i dependence cancels out in _D4, between Eq. (|55l) and 
Eq. ([57)) . These extra pieces therefore do not contain a 
large logarithm and do not need to be exponentiated [see 
the discussion below Eq. (|47| ]. We have again used the 
RPI-III transformation in Eq. ()32p to perform a simple 
check on these expressions. 

E. Electroweak cross section 



as the generations of up and down- type quarks. With 
"channels" we mean the assignment of quark flavor and 
helicity to the lines of the tree-level diagrams in Fig.l. 

The Wilson coefficients Cx are obtained by matching 
at the high scale, running down to the low scale and 
matching at the low scale. 



C{ixi)^D{iii)U{iii,iiH)C{iiH). (59) 



To obtain the electroweak cross section, we need to 
square the matrix elements of the operators in Eq. (1551) , 
sum over quark helicities, flavors, channels and integrate 
over phase space. Schematically, 



1 



2i;2 



XY X2 4 ^ ^ ^ 



hcl. flav. chan. 

X M^i)fA^2) c*xCY{dxY{dY) , (58) 

where 1/(25'^^^^) is the luminosity factor and fi{x) de- 
notes the PDF of flavor i evaluated at momentum frac- 
tion X. The three-body phase space d$3 is given in 
Eq. ([6]), and an identical particle factor of 1/2 (not 
shown) must be included for symmetric phase-space inte- 
grations. With "helicities" we refer to the left- or right- 
handedness of the fermion fields in the operators, which 
runs over {iL, Li?, RK} and fixes the spins of the exter- 
nal particles. Only 1/2 of each of the incoming quarks 
has a given spin, leading to the explicit factor of 1/4. 
The sum over "flavors" runs over {uu, ud, dd] as well 



where U solves the RG equation in Eq. p7|) . In the previ- 
ous sections we described this in detail for the ^-channel 
contribution with incoming quarks, and we will extend it 
to the general case below. To discuss the details, some of 
which are not shown in Eq. (|58p , we consider the helicity 
cases separately. 



1 . Left-Right 



In this case we do not have to worry about identical 
particles. From the i-channel contribution for incoming 
quarks, we can obtain the other channels by a permuta- 
tion of the momenta. For example, if the particles 1 and 
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2 in the original graph were and dj^, we have 



example, 



Ul 


dR 


U , (34)„ , 


Ul 


Ul 


(23),, (243),, 


Ul 


dR 


(24)t , (234)„ , 


dR 


UL 


(12)„,(12)(34)t, 


dR 


UL 


(132)„,(1432)t, 


dR 


dR 


(142),, (1342),, 


Ul 


UL 


(123),, (1243),, 


Ul 


dR 


(13)t,(143)„, 


UL 


dR 


(13)(24)t,(1423)„ 


dR 


UL 


(124)„,(1234)t, 


d-R 


dR 


(14),, (134),, 


dR 


Ul 


(14)(23)t,(1324)„ 



ur , ur : It , (34)„ , 



(60) 



Here we employ the usual notation for elements of the 
permutation group 5*4, e.g. (142) corresponds to 1 ^ 
4 — > 2 — > 1, and 1 denotes the trivial permutation. The 
subscripts on the permutations indicate whether this con- 
tribution corresponds to an s, t or w-channel diagram. 
The permutations are grouped by the type of incoming 
particle, where we (for simplicity) assumed that in the 
original diagram particle 1 and 2 were an up and down- 
type quark, respectively. 

Whether the incoming quarks are up or down- type, 
fixes the flavors of the PDFs except for the generation. 
The electroweak corrections do not depend on the gen- 
eration and the sum over generations differs between the 
various channels. For example. 



Ui,d2 {t,u) : [fu{xi)+fc{xi)][fd{x2)+fs{x2)- 



(61) 



Up to NLL order, we only need the tree-level matrix 
element of the operators at the low scale. For the oper- 
ators in Eq. we find 



{OxViOY) =^SxY X 



ss' for LL, RR 
U'iUji LR, RL 



(62) 



which depends on the helicity of the fields for particle 1 
and 2, as indicated. To obtain the analogue of Eq. (p^ 
for any of the permutations, we may simply permute the 
left and right-hand side of this equation. Note that due 
to the nature of the electroweak corrections, ul and 
are always calculated simultaneously. 



2. Right-Right 

For an up-type and a down-type quark, the approach 
is the same as for LR. If the quarks are of the same 
type, only half the permutation in Eq. (|60l) remain. For 



ur,ur: 



ur,ur: 



ur,ur: 



(23) ,,(234)„, 

(24) « , (243), , 

(13) «,(134),, 

(14) ,,(143)„, 
(13)(24)t,(1324)„. 



(63) 



For identical quarks and helicities, there are interfer- 
ence contributions. Since these are small, we do not in- 
clude them. 



3. Left-Left 

For Oa and Od in Eq. (|33p , for which the flavor types 
are identical, the method for LL is the same as for RR. 
Note that Oc is the (12)(34) permutation of Ob, and Op 
is the (12)(34) permutation of Oe- We may therefore 
use the reduced set of permutations in Eq. (j63l) rather 
than Eq. (pO)) even though for Ob, Oc, Oe and Op the 
flavor types differ. (The reason is that the electroweak 
corrections for ul and dL are always calculated simulta- 
neously.) 



F. Combining electroweak and QCD corrections 

We will now show that we can multiplicatively include 
our electroweak corrections in the QCD cross section. 



O'QCD X ctew 



(64) 



for any given point in phase space and flavor of the in- 
coming partons. Our discussion will be specifically for 
the VBF process but can be fairly straightforwardly ex- 
tended to other cases. Eq. (|64|) is only true up to NLL 
order, since the cross section has the following schematic 
form, 

log a = V i(a,L)" + Lia^Lr + (a,i)" + (a^i)" 



LL NLL 

+ a,(a,L)" -I- ai„(at„L)" + Q!^(q!,L)" + . . 



NNLL 



(65) 



where L denotes a large logarithm and is a weak cou- 
pling (q!2 or ai). At NNLL order, the QCD and elec- 
troweak corrections get mixed and cannot be calculated 
separately. We will see below explicitly where this break- 
down occurs. However, one might only care about getting 
the correct NNLL QCD corrections, since the electroweak 
corrections are smaller. 

To discuss how QCD and electroweak corrections can 
be combined, we will assume that the QCD corrections 
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are also calculated using SCET. The QCD corrections 
only affect the quarks, for which the possible color struc- 
tures are given by 

*37^r^t/(3)*l*47pr^t/(3)*2, *37'^*1*47m*2 . (66) 

In this equation we suppressed the SU{2) structure, 
since the SU{3) and SU{2) generators act in different 
spaces. Including these color structures doubles the ba- 
sis in Eq. (jl4p , which we will view as a tensor product of 
the SU{2) and SU{3) structures. 

At NLL order, the high-scale matching is performed at 
tree-level. Since the diagrams are purely electroweak, it 
trivially factors, 

C{fih) = Csui2) ifJ-h) ® Csu{3) (Ph) ■ (67) 

The one-loop QCD corrections to the high-scale match- 
ing differ between the s, t and u-channel diagrams. Thus 
Eq. ()67|) no longer holds when different channels are com- 
bined, which is inevitable due to interferences between 
channels. 

Up to NLL order, the anomalous dimension completely 
separates 

7 = 7SC/(2) «) lsc/(3) + 1S)7(2) 7SU{3) ■ (68) 

Here Jsu{2) and Jsu(3) denote the usual electroweak and 
QCD anomalous dimensions, i.e. Jsu{2} only acts on the 
SU{2) structures and only involves a2,ai. At one-loop 
order Eq. (|68|) is trivial, because the loop is either QCD or 
electroweak. At two-loop order, Eq. no longer holds 
in general, but it does still hold for the cusp anomalous 
dimension. Since gluons and electroweak gauge bosons 
only couple through fermions and act in different spaces, 
a mixed contribution would require a Cp, with one Cp 
from SU{3) and another from SU{2). However, the two- 
loop cusp does not contain a C|. term. At three-loop 
order, diagrams such as Fig. 2 arise, that produce a C|. 
in the cusp anomalous dimension and violate the form in 
Eq. dMl)- Thus Eq. ^ holds up to NLL order but not 
beyond. Integrating Eq. (|68|) from the high scale to the 
low scale leads to an evolution factor 

u{nh,f^i) ^Usu(2){^^h, ®Usu(3){^^h, , (69) 

which is just the tensor product of the electroweak and 
QCD evolution. 

The low-scale matching only affects the electroweak 
basis, 

^(mO = ■D5c/(2)(/iO «) lsc/(3) ■ (70) 

At one-loop order, -Dsc/(2) is trivially the same as in the 
pure electroweak case. For the large logarithm in the low- 
scale matching it also holds at two-loop order, because 
its coefficient is the cusp anomalous dimension. Since the 
low-scale matching is needed at one lower order than the 
running (see table U), the low scale matching even factors 
at NNLL. 



FIG. 2: Three-loop contribution to the cusp anomalous di- 
mension that violates Eq. (|68[1 . The double lines denote Wil- 
son lines. 

In Sec. Ill El we took the tree- level matrix element of 
the operators at the low scales, which is sufficient for 
the pure electroweak corrections up to NLL order. The 
situation is more complicated for QCD, which becomes 
strongly coupled at low scales. QCD can also have its 
own series of large double logarithms that depend on the 
measurement on the hadronic final state. An example 
of a measurement for an A'^-jet signal at hadron colliders 
is the event shape 7V-jettiness tn [s^. Requiring tat < 
T™* <C 1 defines an exclusive iV-jet measurement and 
induces large QCD logarithms a" log™ tm (with m < 2n) 
from vetoing additional jets. In the presence of QCD 
corrections we still square the operators and get the tree- 
level weights in Eq. (|62|). Since there is no dependence 
on the quark flavor in the perturbative QCD corrections, 
these directly factor. The QCD corrections do depend 
on the quark flavor through the PDFs, so Eq. (|64| only 
holds when the flavor of the incoming quarks is fixed. 
We point out that at one-loop order there are corrections 
involving the gluon PDF, but this contribution is equal 
for all quark flavors and thus preserves the factorization 
of QCD and electroweak corrections. 

III. NUMERICAL RESULTS 

In this section we report our numerical results for the 
integrated cross section. The couplings and parameters 
of the standard model have been set at the electroweak 
scale ^ = Mz according to the data of Ref . [36j . Two- 
loop beta functions are usedfs^ for their running up to 
the high energy scale ^ih- The Higgs mass is assumed to 
be Mh = 125 GeV. 

The low energy matching scale is chosen to be /x; = 
Mz- The cross section should not depend on it, since 
/i; is not a physical scale. We checked that the depen- 
dence is very small because of a cancellation between the 
running and the low-scale matching. However, the elec- 
tromagnetic part of the weak interactions does not get 
integrated out at the low scale, and the cancellation can 
never be perfect. In fact the scale dependence is larger in 
the Left-Right and Right-Right processes where the elec- 
tromagnetic corrections become more important. That is 
not a real issue since the cross section is dominated by the 
Left-Left terms, and the low scale dependence of the to- 
tal cross section is negligible. With our choice ^/ = M^, 
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Ecm [TeV] 

FIG. 3: Tree-level cross section according to Eq. (|58|) but 
with no radiative corrections included, divided by the result 
for Ecm = 2 TeV, as a function of the center-of-mass energy 
Ecm (filled circles). The output of HAWK is reported for 
comparison (open circles). Two sets of cuts have been used: 
PT > 20 GeV, 03 > 10°, 04 < 170° (upper data); pr > 40 
GeV, 6*3 > 30°, 04 < 150° (lower data). 

in the low-scale matching only the coUinear function is 
needed, the remainder being less than 1%. That is be- 
cause the remaining contribution largely cancel. We keep 
all the small terms in our numerical calculation anyway. 
Away from Mz, these terms become important since they 
cancel the low scale dependence of the running. 

As shown in table U at NLL order only the tree 
level high energy matching is required, and the one-loop 
matching terms can be neglected. As for any tree-level 
calculation, the result does depend on the high energy 
scale fj,h because there is no cancellation between the run- 
ning and the high energy matching at tree-level. Thus, 
as usual, the high scale must be thoroughly chosen in or- 
der to keep the neglected one-loop terms small enough. 
Typically, the one-loop matching terms contain the log- 
arithmic terms lnsij//i^, thus whenever the Sij are of 
the same size, we take /i^ « Sij and the one-loop match- 
ing terms can be safely neglected, as they do not contain 
large logarithms. In VBF we are mainly interested in the 
t channel, and the relevant logarithms that appear in the 
one-loop matching are ln(— taZ/i^) and ln(— ^4//!^). How- 
ever there are small phase space regions where ^3 ^ 
or t4 ^ ^3 , and the choice of the correct high scale turns 
out to be ambiguous. The best way out would be a full 
calculation of the one-loop high energy matching terms, 
but that would be beyond the aim of the present pa- 
per. Moreover, avoiding difficult electroweak loop calcu- 
lations was one of the motivations of the present study. 
It is worth mentioning that for VBF the same ambigu- 
ous choice of the high scale is also encountered in the tree 
level calculation of the cross section whenever ^3 and t4 
are not of the same size. 

As a good interpolation between phase space regions 
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FIG. 4: Single terms contibuting to the tree-level cross section 
as a function of the center-of-mass energy -Ecm. Cuts are: 
PT > 20 GeV, 6I3 > 10°, 04 < 170°. 



where ^3 sa and regions where one of the variables 
is small, we set the high scale to the minimum between 
Mz and the geometric average fif^ = \/t^tl. At this scale 
the sum of logarithmic terms ln(— ^3//^^) + ln(— t4//i^) 
reaches its minimum in the one-loop matching. More- 
over this choice has the merit of stopping the running 
whenever one of the Mandelstam variables is too small. 
We will discuss the sensitivity of the result to the choice 
of the high scale, and show that it is comparable to the 
sensitivity of the standard tree-level cross section. ^ 

The aim of the paper is a test of the method by a com- 
parison with the one-loop fixed order perturbative calcu- 
lation of the code HAWK that is based on the work of 
Ref . [H, [l^ . In our comparison we adopt larger cuts than 
phenomenologically required, in order to be sure that we 
are on safe grounds for the SCET approximation. We 
use CTEQ6 PDFsji^l ^"^^ neglect the very small con- 
tribution of t and b quarks. Moreover in both codes we 
neglect the s-channel contribution and interference terms 
that are known to be small with VBF cuts. 

At tree level we find a fair agreement between the 
codes, provided that the masses of the vector gauge 
bosons are restored in the denominators of the coefficients 
in Eq. (|16p . While that is not so relevant when relative 
cross sections are reported with large cuts, the effect of 
masses becomes crucial when the total cross section is 



^ We would like to stress once more that, at variance with the vari- 
ational approaches to the electroweak sector [38h4C| . the principle 
of minimum sensitivity here does not apply even if /i/j is not a 
physical scale. In fact the sensitivity only measures the depen- 
dence on /i/i of the neglected terms, and this dependence might 
be high even when the neglected terms are small. The best choice 
of /ift should make the omitted terms small, as for the method 
of minimal variance l4lll . 
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FIG. 5: Relative electroweak correction [(Jew —o'TREE\/(^tree 
as a function of the center-of-mass energy i?cm evaluated by 
Eq. (|58|l in this paper (filled circles), compared with the 
output of HAWK (open circles). Cuts are: pr > 20 GeV, 
e-i > 10°, 6I4 < 170°. 
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FIG. 6: The K-factor K — uew /lytree as a function of 6 — 
03 = 180° - 04, at = 90°, XI = X2 ^ 0.5, Ecu, = 14 TeV 
and E3 = xiX'zEcm/^. 



evaluated by integrating over the whole phase space. In 
relative terms, the tree-level cross section, without any 
radiative correction, is reported in Fig. 3 for several values 
of the center-of-mass energy Ecm and for different cuts, 
as a test of the integration over the PDFs according to 
Eq. dSll). 

Hereafter we adopt the following choice of cuts on 
angles and transverse moment: 6*3 > 10°, 64 < 170°, 
Pt > 20 GeV. The single terms contributing to the cross 
section at tree-level are reported in Fig. 4. The ud -> du 
process is the dominating one, followed by the other Left- 
Left terms. The Left-Right terms are two order of magni- 
tude smaller, while the Right-Right terms are negligible 
and have not been included in the calculation. A sum 
over generations of the up- and down-type quarks is in- 
cluded through the PDFs in the terms of Fig. 4. 

The electroweak corrections have been evaluated inte- 
grating in Eq. (j58p with the same cuts. In order to com- 
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7: The K-factor K = uew /o'tree as a function of ip at 
10°, 64 = 170°, xi = X2 = 0.5, Ecm = 14 TeV and 

XlX2Ecin/2. 
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pare with HAWK no QCD corrections have been included 
in both codes, and only virtual electroweak correction are 
considered. We find large electroweak corrections in wide 
sectors of the phase space, expecially for the Left-Right 
processes that get suppressed by even more than 90% 
at some spots (see Figs. 6 and 7). However, after inte- 
grating with PDFs, the corrections are not dramatic. In 
fact the corrections are smaller for the ud — du process 
which dominates in the cross section, while the highly 
suppressed terms, like ULdn, have a very small weight on 
the total result. 

In Fig. 5 we show the relative electroweak correction, 
defined as the ratio 



a 



^tree 



(71) 



as a function of the center-of-mass energy. The correc- 
tion is less than 3% at Ecm = 2 TeV, in good agreement 
with the output of HAWK. At larger energies the cor- 
rection grows faster than predicted by HAWK, reaching 
9% at the LHC energy Ecm — 14 TeV, to be compared 
with 5% predicted by HAWK. At the higher energies, the 
larger contribution to the electroweak correction comes 
from the running of the coefficients by the anomalous 
dimension. Thus we argue that fixed-order perturbative 
calculations might miss part of the correction at the LHC 
energy scale. 

It is instructive to see how the single terms behave in 
the phase space, and for each of them we define a sort of 
K-factor as 



K 



(7tree{0Z,S4,(p) ' 



(72) 



where the differential cross sections are evaluated for 
fixed values of xi,X2, E^ and Ecm, and for a given set 
of angles. In the physically relevant range of the param- 
eters we do not find any important dependence on the 
angles, and the K-factor is almost flat at small angles. 
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FIG. 8: The relative change of the cross section in the process 
ud — >■ du is reported for a 10% increase of the high scale fih 
(filled circles). The tree- level result is shown for comparison 
(open circles). Cuts are the same as in Fig. 5. 



parison, and is evaluated by taking fj,h as renormalization 
scale. The tree-level cross section depends on the high 
scale because of the running of the couplings, and we 
find a steady -0.6% change for all the energies. As shown 
in Fig. 8, when the electroweak corrections are included 
the cross section dependence increases with the increase 
of energy, going from -1.2% to -1.8%. This dependence 
is also negative, and remains comparable in size to the 
tree-level value. Moreover these dependences cancel each 
other in the relative electroweak correction Aa/a which 
we reported in Fig. 5. In fact, differenziating Eq. ((7T|) and 
assuming atree ~ ctew 



a 



Sctew Satre 



(73) 



Thus the overall effect on the relative correction goes 
from -0.6% to -1.2% which seems to be more than ac- 
ceptable if compared with the natural -0.6% change of 
the tree-level cross section. 



IV. CONCLUSIONS 



For instance in Fig. 6 we show the dependence on 03, 6^. 
The K-factor is reported as a function of 6* = 6'3 with 
04 = TT — 03 and = 7r/4 at a typical set of parameters: 
Ec-m — 14 TeV, xi = X2 = 0.5 and = a;ia;2£'cm/2. 
For 9 < 90° the dependence on 9 is very small and is not 
expected to have any important relevance on the angular 
distribution of the jets. It is remarkable the very large 
suppression of the ulcIe. process which has no relevant 
effect on the total cross section that is dominated by the 
Left-Left processes. The dependence on the azimuthal 
angle ip seems to be even smaller, but becomes more rel- 
evant when ip approaches 180°, indicating that at large 
values of (y9 w 180° the cross section gets enhanced by the 
electroweak corrections, as predicted by the NLO calcu- 
lation of Ref.[19|. For instance in Fig. 7 the dependence 
on ip is shown at small angles 03 = 10°, 04 = 170° and 
for the same set of parameters as before. 

As discussed above, our result must depend on the high 
scale /ift, because the tree-level matching cannot cancel 
the dependence. The same problem arises in the calcu- 
lation of the tree-level cross section, and one takes for 
granted that the choice of the high scale must keep the 
logarithms small in the higher order terms of the per- 
turbative expansion. However, even if our choice for the 
high scale seems to work well, it would be desirable to see 
how the results depend on the choice of fifi, expecially at 
low energies where the Mandelstam variables might be 
comparable with the low scale = Mz- We can check 
this dependence on the most dominant term in the cross 
section, namely the ud — > du process, and we can com- 
pare with the dependence of the tree-level cross section. 

In Fig. 8 the relative change 5a /a is reported for a ten 
per cent increase of the high scale S^h/lJ-h = 0.1. The 
change of the tree-level cross section is reported for com- 



Even if VBF is the second largest production channel 
for the Higgs boson, it is a clean and pure electroweak 
process, and its precise measurement is very important 
for constraining the Higgs couplings and thus identifying 
the nature of the Higgs sector. A full control of the ra- 
diative corrections might be also relevant for ruling out 
extensions of the standard model or exo tic p rocesses that 
have been predicted in the Higgs sector [43l li^l 

We have shown that SCET and the method of Refs. [1, 
20] may be used to directly obtain electroweak correc- 
tions for generic processes without the need to perform 
difficult electroweak loop calculations that would have 
to be carried out for each individual process. The elec- 
troweak corrections have been obtained in analytical 
form, and can be easily inserted in the software pack- 
ages that have been developed for computing QCD cor- 
rections to cross sections. This is one of the first tests 
of the formalism at the level of a hadronic cross section, 
and demonstrates the viability of the method that has 
the merit of providing a resummation of the large dom- 
inating logarithms at LHC energies. Moreover, in the 
special case of VBF, the amplitude is explicitly propor- 
tional to the VEV, so the effective field theory operator is 
not a gauge singlet and standard resummation methods 
do not apply. Thus VBF provides a very interesting test 
of the method, besides the phenomenological relevance 
of the process in the study of the Higgs boson at LHC. 

In our test we used the extension to the VBF process 
as derived in Ref . [l2| , and compared the integrated cross 
section with the NLO output of the code HAWK, based 
on the calculation of Ref. [isl . [Toj . 

We find large electroweak corrections in wide sectors of 
the phase space, expecially for the Left-Right processes 
where the corrections reach 90%. They mostly arise from 



15 



the RG running of the Wilson coefficients of the effec- 
tive operators. However, after integrating with PDFs, 
the corrections are not dramatic. In fact the corrections 
are smaller for the ud — ?> du process which dominates in 
the cross section, while the highly suppressed terms, like 
ULdR, have a very small weight on the total result. 

Below 2 TeV our results are in fair agreement with 
HAWK, while at higher energies the SCET formalism 
predicts corrections that are slightly larger, and grow 
up with increasing the center-of-mass energy. Since the 
larger effect comes from the running, we argue that fixed- 
order perturbative calculations might miss part of the 
correction at the LHC high energy scale. We also ex- 
plored the angular dependence of the electroweak correc- 
tion which is small and in general agreement with NLO 
calculations [H [Til. 

While we checked that the sensitivity to the high scale 
/i/i is not too high, and is comparable to the sensitivity of 
the tree-level cross section, we are aware that the correc- 
tions must depend on the choice of /i/i because no cancel- 
lation is provided by the tree-level high scale matching, 
and that a thoroughly choice for is important. Our 
choice of a geometric average /i^ = V^sti seems to work 
well at low energies, and is expected to work even bet- 
ter at higher energies since the Mandelstam variables get 
much larger than the electrowek scale Mz- A full one- 
loop calculation of the high energy matching would be 
required in order to test the choice, but that goes out 
of the aim of the present paper and will require further 
work. 



given by 

Ao{m) 

ABq^ui, s) 



m" 1 — loa 



TO'' 

log — 

4z 



2V4z- Itan"^ 
1 



1 



; tan 



-II, (A2) 



where z — m} js. These functions are real if s < 4to^, 
while for s > Am? they acquire an imaginary part and 
can be explicitly written as 



BQ{Tn,s) = 2-log 



ABq^iti, s) 



TO 

2z 



, l + ^/T^ 

log , —ITT 



l-^/T^ 



log ^= - 1 



27ri z 

sVl - 4z ■ 
(A3) 



Appendix B: Explicit form of Higgs wave function 
R'' and tadpole r'' 

The finite part of the Higgs wave function i?'* and the 
scalar tadpole graph F'* appear in the low-scale matching 
in Eqs. (|48| and ([54|. For completeness we give their 
explicit expressions here, which can for instance be found 
in Ref. [iq. 

In gauge, the Higgs wave function 6R'^ = i?'* — 1 
is given by 
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Appendix A: Passarino-Veltman functions 

The Passarino-Veltman equal-mass functions [i^ Aq, 
Bq and ABq are given by the following analytical integral 
representations (with the normalization of Ref. (46t ) 

f d^fc 1 
167r2 ~ J (2^)'* fc2 _ ^2 + iO ' 

Ba[m,p ) 



IGtt 



{2ttY (fc2-TO2+iO)[(fc+p)2-TO2+iO] 



d 



ABa{m,s) = —Ba{m,s) 



(Al) 



In the MS scheme, the finite part of these functions are 



1 



167r2ti2 



2Af|Bo(Mz, Ml) + AM^BoiMw, Ml) 
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MtAB^{Mh, Ml) - 2m^MmuMl) 



2MlMl 



1 



Mt-&M^z]ABo{Mz,Ml) 



{AMlM^ - M^ ~ 12M^)ABo{Mw,Ml) 
{2Mlm',-8mt)ABo{muMl) 



(Bl) 



Here, nit is the mass of the top quark, and the contribu- 
tion from the other fermions have been neglected. The 
Passarino-Veltman functions Bq and ABq are given in 
App. 13 Note that while only the real part of the wave 
function can enter in the renormalization, the full com- 
plex value must be kept in the low-scale matching. 

The tadpole contribution F'' is related to the shift Sv of 
the VEV. By a self-consistent evaluation it can be written 
in gauge as 



Mlv 



Sv 

V 



1 



167r2i;2M2 



2AC 



4M^ 



(B2) 



^MlAoiMh) - (3M| + ^a4)Ao{Mz) 



(6M^ + MI)Aq{Mw) + AmlAo{mt) 
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where again only the fermionic term for the top quark has been retained. 
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